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Abstract 

All the segments (or their continuations) of a billiard trajectory inside an ellipsoid of R" are 
tangent to n — 1 quadrics of the pencil of confocal quadrics determined by the ellipsoid. The 
quadrics associated to periodic billiard trajectories verify certain algebraic conditions. Cayley 
found them in the planar case. Dragovic and Radnovic generalized them to any dimension. 
We rewrite the original matrix formulation of these generalized Cayley conditions as a simpler 
polynomial one. We find several remarkable algebraic relations between caustic parameters and 
ellipsoidal parameters that give rise to nonsingular periodic trajectories. 
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1. Introduction 

One of the best known discrete integrable system is the billiard inside ellipsoids. All the 
segments (or their continuations) of a billiard trajectory inside an ellipsoid of R" are tangent 
to n — 1 quadrics of the pencil of confocal quadrics determined by the ellipsoid [1-3]. This 
situation is fairly exceptional. Quadrics are the only smooth hypersurfaces of R", n > 3, that have 
caustics flm. A caustic is a smooth hypersurface with the property that a billiard trajectory, 
once tangent to it, stays tangent after every reflection. Caustics are a geometric manifestation of 
the integrability of billiards inside ellipsoids. 

Periodic trajectories are the most distinctive trajectories, so their study is the first task. There 
exist two remarkable results about periodic billiard trajectories inside ellipsoids: the generalized 
Poncelet theorem and the generalized Cayley conditions. 

A classical geometric theorem of Poncelet JaLZI] implies that if a billiard trajectory inside an 
ellipse is periodic, then all the trajectories sharing its caustic are also periodic. Its generalization 
to the spatial case was p roved by Darboux Jfl]. The extension of this result to arbitrary dimensions 
can be found in 19 H12I1 . The generalized Poncelet theorem can be stated as follows. If a billiard 
trajectory inside an ellipsoid is closed after mo bounces and has length Lq, then all trajectories 
sharing the same caustics are also closed after mo bounces and have length Lq. Thus, a natural 
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question arises. What caustics do give rise to periodic trajectories? The planar case was solved 
by Cayley lfl3i[l4ll in the XIX century. Dragovic and Radnovic Il5l[l6ll found some generalized 
Cayley conditions for billiards inside ellipsoids fifteen years ago. They have also stated similar 



conditions in other billiard frameworks; see [ 17-21]. 



For simplicity, let us focus on the spatial case. Let Q : x /a + y /b + z /c — 1 be the triaxial 
ellipsoid with ellipsoidal parameters < c < b < a. Any billiard trajectory inside Q has as 
caustics two elements of the family of confocal quadrics 

( x 2 V 2 7 2 

Qa - \ (x, y, z) e R 3 : + ■ + = 1 

[ a — A b — A c — A 

We restrict our attention to nonsingular trajectories. That is, trajectories with two different caus- 
tics which are ellipsoids: < A < c, hyperboloids of one sheet: c < A < b, or hyperboloids 
of two sheets: b < A < a. The singular values A e {a, b, c} are discarded. There exist some 
restrictions on the caustics and Q^. It is known that EH1, H1H1, EH2, and H1H2 are the 



only feasible caustic types; see 1221 123H . The meaning of these notations is evident. 

The generalized Cayley condition in this context can be expressed as follows. The billiard 
trajectories inside the triaxial ellipsoid Q sharing the caustics and Qx 2 are periodic with 
elliptic period m > 3 if and only if 

( /m + l ■•• U ^ 



rank 



flm—X ' ' ' fm+2 



< m — 2, 



where /(f) = £/> /z*' : = VP - t/a)(l - t/b)(l - f/c)(l - fMiXl - t/A 2 ). We claim that most 
results related to generalized Cayley conditions become simpler when expressed in terms of the 
inverse quantities I /a, l/b, l/c, l/A\, and l/^- Proposition[T2lis a paradigmatic example. 

The elliptic period is defined in Section [2] Roughly speaking, the difference between the 
period mo and the elliptic period m of the periodic billiard trajectories sharing two given caustics 
is that all those trajectories close in Cartesian (respectively, elliptic) coordinates after exactly mo 
(respectively, m) bounces. We will see that either m = mo/2 or m = mo. 

The previous matrix formulation is very nice from a theoretical point of view, but it has strong 
limitations from a computational point of view. We will see in Section [3] that it can be written 
as a system of two homogeneous symmetric polynomial equations with rational coefficients of 
degrees m 2 - 2 and m 2 - 1 in the variables l/a, l/b, l/c, 1/A\, and I//I2. Thus, both degrees 
grow quadratically with the elliptic period m, which turns this approach into a tough challenge. 
In particular, to our knowledge, the caustic parameters A\ and A% have never been explicitly 
expressed in terms of the ellipsoidal parameters a, b, and c for any m > 3. 

We will rewrite this matrix formulation as a computationally more appealing one which gives 
rise to (non-symmetric) homogeneous polynomial equations whose degrees are smaller than the 
elliptic period m. We will find the following remarkable algebraic relations between caustic and 
ellipsoidal parameters using the new formulation. The billiard trajectories inside the ellipsoid Q 
sharing the caustics and Qx, are periodic with: 

• Elliptic period m = 3 if the roots of f 3 - (f - c)(t - b)(t - a) are the caustic parameters; 

• Elliptic period m = 4 if there exists d e R such that the roots of f 4 - (f - d) 2 (t - Ai)(t - A2) 
are the ellipsoidal parameters; and 
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• Elliptic period m — 5 if the roots of t 5 - (t - c)(t - b)(t — a)(t - Ai)(t - Az) are double ones. 

Let us compare both formulations in the third case. On the one hand, the two homogeneous 
symmetric polynomial equations obtained from the matrix formulation have degrees 23 and 24 in 
the variables l/a, l/b, 1/c, l/Au and l/A 2 . On the other hand, t 5 -(t-c)(t-b)(t-a)(t-Ai)(t-A 2 ) 
is a polynomial of degree four in a single variable. Clearly, the polynomial formulation leads to a 
much simpler problem. Nevertheless, it should be stressed that the matrix formulation determines 
all periodic billiard trajectories with elliptic period m = 5. On the contrary, we find just a subset 
of such trajectories using the polynomial f 5 - (f - c)(f - b)(t - a){t - A\){t - Aq). 

Another natural question about periodic billiard trajectories is the following one. Which are 
the triaxial ellipsoids of R 3 that display periodic billiard trajectories with a fixed caustic type and 
a fixed (elliptic) period? A numerical approach to that question was considered in [24], where 
the authors computed several bifurcations in the space of ellipsoidal parameters. We will find the 
algebraic relations that define the bifurcations associated to small elliptic periods. For instance, 
we will see that there exist periodic billiard trajectories with elliptic period m — 3 and caustic 
type EH1 if and only if c < ab/(a + b + ^fab). 

For brevity, we will not depict billiard trajectories inside triaxial ellipsoids of R 3 . The reader 
interested in 3D graphical visualizations is referred to [25], where several periodic billiard tra- 
jectories with small periods are displayed from different perspectives. 

We complete this introduction with a note on the organization of the article. In Section [2] 
we review briefly some well-known results about billiards inside ellipsoids, recalling the matrix 
formulation of the generalized Cayley conditions obtained by Dragovic and Radnovic. We also 
introduce the concept of elliptic period. The practical limitations of the matrix formulation are 
exposed in Section [3] Next, we present the polynomial formulation in Section [4] In Section [5] 
we carry out a detailed analysis for minimal elliptic periods, whereas the study of more general 
elliptic periods is postponed to Section [6] The previous results are adapted to billiards inside 
ellipses of R 2 and inside triaxial ellipsoids of M 3 in sections Q and [8] respectively. 

2. Preliminaries 

In this section we recall several classical results and their modern generalizations about bil- 
liards inside ellipsoids that go back to Jacobi, Chasles, Poncelet, Darboux, and Cayley. 
We consider the billiard dynamics inside the ellipsoid 



The degenerate cases in which the ellipsoid has some symmetry of revolution are not considered 
here. This ellipsoid is an element of the family of confocal quadrics 



We note that Q,\ = for A > a„. Thus, there are exactly n different geometric types of nonsingular 
quadrics in this family, which correspond to the cases 




< a i < • ■ ■ < a n . 



(1) 




A €(-00,0!), /le(ai, 02), Ae(a n -i,a„). 
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For instance, the confocal quadric Q\ is an ellipsoid if and only if A e (-oo,^). On the other 
hand, the meaning of Q\ in the singular cases A e {a\, . . . , a n ) is 



Q a . =Hj = [x=(x u ,.., x n ) e R" : Xj = o) . 
The following theorems of Jacobi and Chasles can be found in HJ-dl . 



Theorem 1 (Jacobi). Any generic point ieM" belongs to exactly n distinct nonsingular quadrics 
Quo, ■ ■ ■ , Qfin-i such that Ho < a\ < /ii < a 2 < ■ ■ ■ < a n -\ < p n -\ < a n . 

We denote by fj. = Qio, , . . ,ju„_i) e M." the Jacobi elliptic coordinates of the point x = 
(*!,..., x^. Cartesian and elliptic coordinates are linked by relations 

2 nEo'^-A*,-) . , 

xn = , j = l,...,n. 

Hence, a point has the same elliptic coordinates that its orthogonal reflections with respect to the 
coordinate subspaces of M". A point is generic, in the sense of Theorem Q] if and only if it is 
outside all coordinate hyperplanes. When a point tends to the coordinate hyperplane Hj, some 
of its elliptic coordinates tends to aj. 

Theorem 2 (Chasles). Any line in W is tangent to exactly n—l confocal quadrics , . . . , Qa,,^ ■ 

It is known that if two lines obey the reflection law at a point x e Q, then both lines are 
tangent to the same confocal quadrics. Thus, all lines of a billiard trajectory inside the ellipsoid 
Q are tangent to exactly n—l confocal quadrics Qa„-, , which are called caustics of the 

trajectory, whereas A\, . . . , A n -\ are the caustic parameters of the trajectory. We will say that a 
billiard trajectory inside Q is nonsingular when it has n—l distinct nonsingular caustics. We 
mostly deal with nonsingular billiard trajectories in this paper. 

The caustic parameters cannot take arbitrary values. For instance, a line cannot be tangent 
to two different confocal ellipsoids, and all caustics parameters must be positive. The following 



complete characterization was given in [22, 23]. 

Proposition 3. Let A\ < ■ ■ ■ < A n -\ be some real numbers such that 

{fli,. ..,«„} n [A\,...,A n -\) = 0. 

Set flo = 0. Then there exist nonsingular billiard trajectories inside the ellipsoid Q sharing the 
caustics Qa x , ■ . ■ , Qa„^ if and only if 

A k € (a k -i,a k ) U (a k , a k+1 ), k=\,...,n-\. (2) 

Definition 1. The caustic type of a nonsingular trajectory is the vector g — (q\, . . . , g n -\) £ 
such that 

A k e (a n , a ft +i), k = 1, . . . ,n — 1. 

We know from Proposition[3]that g k e {k - 1 , k] for k — 1, . . . , n — 1. Hence, there are exactly 
2"~' different caustic types. The two caustic types in the planar case correspond to ellipses: 
g\ = 0, and hyperbolas: g\ = 1. The four caustic types in the spatial case correspond to EH1: 
g = (0, 1),H1H1: g = (1,1), EH2: g = (0,2), andHlH2: g= (1,2). The notations EH1, H1H1, 
EH2, and H1H2 were described in the introduction. 

Next, we recall a result about periodic billiard trajectories inside ellipsoids. 
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Theorem 4 (Generalized Poncelet theorem). If a nonsingular billiard trajectory is closed after 
mo bounces and has length Lq, then all trajectories sharing the same caustics are also closed 
after mo bounces and have length Lq. 

Poncelet proved this theorem for conies |@]. Darboux generalized it to triaxial ellipsoids of 
R 3 ; see |@]. Later on, this result was generalized to any dimension in 1I9I- 121. 



The periodic billiard trajectories sharing the same caustics also have the same winding num- 
bers. In order to introduce these numbers, we set 

{ci,...,c 2n -i} = {a\,...,a n } U [A u . ..,A n -i}, (3) 



and f{t) = -^n?=i '(1 _ t/cj). We deal with nonsingular billiard trajectories inside ellipsoids 
without symmetries of revolution, so the parameters c\, . . . , c 2n -\ are pairwise distinct, and we 
can assume that cq := < c\ < ■•■ < c 2n -\. 

Theorem 5 (Winding numbers). The nonsingular billiard trajectories inside the ellipsoid Q 
sharing the caustics Q\ x , . . . , Q\ n _ x are periodic with period mo if and only if there exist some 
positive integer numbers mi, ... , m„_i such that 

n-l 

(4) 



2J-\) } mj d? = 0, Vz = 0,...,n 

j=o J°v 



Each of these periodic billiard trajectories has mj points at Q Clj and mj points at Q C2 m ■ Besides, 
{c2j, C2j+i} H {a\, . . . , a n } => mj even. Finally, gcd(mo, . . . , m„_i) € {1,2}. 

Let Lq be the common length of these periodic billiard trajectories. Let x(i) be an arc-length 
parametrization of any of these trajectories. Let p(t) — (po(t), . . . ,//„_i(f)) be the corresponding 
parametrization in elliptic coordinates. Then: 

1. c 2 j < pj(t) < c 2 j+\ for all t e R 

2. Functions pj(t) are smooth everywhere, except /io(f) which is non-smooth at impact points 
— that is, when :c(f*) e Q — , in which case /ig(f*+) = — — ) + 0. 

3. If p. jit) is smooth at t — t+, then l^'j(t*) = 0o Pj(t*) e [c 2 j, Qj+i). 

4. Pj(f) makes exactly mj complete oscillations (round trips) inside the interval \_Czj, C27+1] 
along one period < t < Lq. 

5. p(t) has period L — LqI gcd(mo, . . . , w„_i). 

Definition 2. The numbers mo, . . . , m n -\ are called winding numbers. Theorem|5]contains three 
equivalent definitions for them: by means of the property regarding hyperelliptic integrals given 
in (01), as a geometric description of how the periodic billiard trajectories fold in R", and as the 
number of oscillations of the elliptic coordinates along one period. 

Most of the statements of Theorem [5] can be found in iTTsL I20I1 . but the one about the even 
character of some winding numbers and the ones regarding gcd(mo, . . . ,m„_i). The first state- 
ment is trivial; suffice it to realize that a periodic billiard trajectory can only have an even number 
of crossings with any coordinate hyperplane. The second ones follow from the oscillating be- 
haviour of elliptic coordinates along billiard trajectories described in Theorem [5J suffice it to 
note that all elliptic coordinates make an integer number of complete oscillations inside their 
corresponding intervals along one half-period Lq/2 when gcd(mo, . . . , m„_i) = 2. 



The following conjecture was stated in [24], where it was numerically tested. 
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Conjecture 1. Winding numbers are always ordered in a strict decreasing way; that is, 



2 < m n -\ < ■ ■ ■ < ni\ < mo- 



lt is known that the conjecture holds in the planar case. If this conjecture holds, then any 
nonsingular periodic billiard trajectory inside Q has period at least n + 1. By the way, there are 
periodic billiard trajectories of smaller periods, but all of them are singular — they are contained 
in some coordinate hyperplane or in some ruled quadric of the confocal family. 

In light of the last item of Theorem[5] we present the following definitions. 



Definition 3. The elliptic period m and the elliptic winding numbers mo, . . . 
gular periodic billiard trajectories with period mo and winding numbers mo, 



m„-\ of the nonsin- 
. . , m„_i are 



m — mo/d, mj—mj/d, 



where d = gcd(mo, . . . , m„_i). 



Roughly speaking, the difference between the period mo and the elliptic period m is that pe- 
riodic billiard trajectories close in Cartesian (respectively, elliptic) coordinates after exactly mo 
(respectively, m) bounces. In order to clarify this difference, let us consider the six planar peri- 
odic trajectories shown in Figure[TJ see Section|7] Only the trajectory in Figure [T(c")| verifies that 
m = mo- On the contrary, the trajectories in figures [T(a)]|l(b)| and | 1 (e)| (respectively, Figure [T(d)] i 
(respectively, Figure [T(T)] i have even period mo and any of their impact points becomes its reflec- 
tion with respect to the origin (respectively, the vertical axis) (respectively, the horizontal axis) 
after mo/2 bounces, so they have elliptic period m = mo/2. 

It turns out that given any ellipsoid of the form ([TJ and any proper coordinate subspace of 
M", there exist infinitely many sets of n — 1 distinct nonsingular caustics such that their tangent 
trajectories are periodic with even period, say mo, and any of their impact points becomes its 
reflection with respect to that coordinate subspace after mo/2 bounces. We will not prove this 
claim, since the proof requires some convoluted ideas developed in 12-4 12511 . 

It is natural to look for caustics giving rise to periodic billiard trajectories inside that ellipsoid. 
Such caustics can be found by means of certain algebraic conditions, called generalized Cayley 
conditions. They are found by working in elliptic coordinates, so they depend on the elliptic 
period m, not on the (Cartesian) period mo. 



Theorem 6 (Generalized Cayley conditions). The nonsingular billiard trajectories inside the 
ellipsoid Q sharing the caustics ,...,Qx 
m > n and 

fm+\ 



rank 



are periodic with elliptic period m if and only if 

fn+l 



f: 



2m- 1 



< m — n + 1, 



m+n- 1 / 



where f(t) = £/>o//f' := ^niV'd - t/a). 



Cayley proved this theorem for conie s 1 1 3 f l - Later on, this result was generalized to any 
dimension by Dragovic and Radnovic in jl5ill6ll . These authors have also given similar Cayley 
conditions in many other billiard frameworks; see Il7l42in . 



Definition 4. C(m, n) denotes the generalized Cayley condition that characterizes billiard trajec- 
tories of elliptic period m inside ellipsoids of R" given in Theorem[6] 
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3. On the matrix formulation of the generalized Cayley conditions 



The matrix formulation of the generalized Cayley condition C(m,n) stated in Theorem|6]is 
very nice from a theoretical point of view, but has strong limitations from a practical point of 
view. Let us describe them. 

The function f(f) is symmetric in the inverse quantities y,- = 1/c,-. In order to exploit it, we 
introduce some notations about symmetric polynomials. Let Q*} om ' sym [x\, . . . , x s ] be the vectorial 
space over Q of all homogeneous symmetric polynomials with rational coefficients of degree / 
in the variables x\, . . . , x s . Let &i(xi , . . . , x s ) be the elementary symmetric polynomial of degree / 
in the variables X\, . . . , x s . That is, n/=i(l + x i0 = 2/>o e /(*i> ■ • ■ > x s )r, so e/(xi, ...,x s ) = for 
all I > s. Clearly, e/ = &i{x\, . . . , x s ) e (Q)! lom ' sym [xi , . . .,x s ]. 

We stress that // = fi{y\, . . . , yin-i) £ Q?° m,sym [')'i, ■ • ■ , yin-i], which is one of the reasons for 
the introduction of the inverse quantities y,- = 1/c,. Indeed, using that f 2 (t) = ]~[f=i f'(l _ 7/0> we 
get the recursive relations 

i-i 

fo = l, 2/, = (-l)W,-.-,72 n -i) -£/*//-*, W>1. 

Hence, it is possible to compute recursively all Taylor coefficients //, although their expressions 
are rather complicated when / is big. Nevertheless, the computation of the Taylor coefficients 
fn+\, ■ ■ ■ ,fzm-\ is me simplest step in the practical implementation of the generalized Cayley 
condition C(m, ri). Next, we must impose that all (m - n + 1) x (m - n + 1) minors of the matrix 
that appear in Theorem [6] vanish. For simplicity, let us consider the minors formed by the first 
m - n rows and the (m -n + l)-th row of that matrix, for I = 1 , . . . , n— 1 . Then the Cayley condition 
C(m, n) can be written as the system of n - 1 polynomial equations 



Mm,n,l = Mmjnjfyx,. . . , yin-l) - 



fm+l 

flm-n 
flm-n+l 



fn+l 

fn 
fm+l 



0, !</<«-!. 



(5) 



One can check that M m , nJ (y u ..., y 2 „-\) e Q^+^ m _ n+ ,[yi, . . . ,72»-i] from the Leibniz formula 
for determinants. This implies that the resolution of system © is a formidable challenge, even 
from a purely numerical point of view and for relatively small values of m. 

We want to write down the solutions of system © in an explicit algebraic way. Let us focus 
on the planar case n = 2, when condition C(m, 2) becomes a single homogeneous symmetric 
polynomial equation of degree m 2 - 1 in three unknowns; namely, 



M m = M m (yi,y 2 ,y 3 ) := 



m+l 



flm-l 



h 

fm+l 



For instance, condition C(2, 2) can be easily solved, since 



-16M 2 = y\ + y\+y\- y 2 ^ - y 2 y 3 - y\y\ - yjy 3 - y\y\ - y\y 2 + lyiyiy^ 

= (71 - 72 - 73)(73 - 71 - 72X72 - 73 - yi)- 
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The inverse quantities y, = 1 /c, verify that < y$ < y% < y\, since < c\ < c 2 < C3. Therefore, 
only the first factor of the above formula provides a feasible solution, and so condition C(2, 2) 
has a unique solution: 

7i = 72 +73- 

The computations for condition C(3, 2) are much harder, so we have implemented them using 
a computer algebra system. We got the factorization -16384M3 = qoq\q 2 qi, where 

<7o = j\ +72 + 73 - 27172 -27173 -2y 2 y 3 , 

q k = lyl-liyi + y^yu-iyi-yjj 1 , (i, j, k} = {1,2,3}. 

The first factor go can, in its turn, be factored as 

qo = ( VtT - V72 - VtIX VtT + V72 - VtIX VtT - VtT + VtIX VtT + VtT + Vrl)- 

The factor qo provides a unique feasible solution: yjyi = yjyi + yjyi, because < 73 < 72 < 7i ■ 
Next, we consider the factor qk as a second-order polynomial in the variable % with coefficients 
in Z sym [yi,yj]. Then we get the solutions 

±, \ 7; + 7; 2 r~ 2 : 

7k = 7 k (7h 7j) ■= — — ± 3 A/rf + 7J - 7i7j- 

It turns out that y k < < max(y,, yj) < y^ , so only the factor qi gives a solution compatible with 
the ordering < 73 < 72 < 7\\ namely, y\ - y\{yi, 73). Hence, C(3, 2) has only two solutions: 



VrT = Vr? + Vrl and 3yi = y 2 + 73 + 2 Jy^ + y\ - y 2 y3- (6) 

We have tried to write down explicitly the solutions of system (|5]l in other cases, but we did 
not succed, even after implementing the computations in a computer algebra system. This shows 
the limitations of the matrix formulation. 



4. A polynomial formulation of the generalized Cayley conditions 

Theorem 7. Let r(i) — n?=i '(1 _ f / c /') and f{t) — y/r(f). The generalized Cayley condition 
C(m, n) is equivalent to any of the following two conditions: 



(i) There exists a non-zero polynomial s(t) G R m _„[f] such that 

{s(t)f(t)} = 0, l = m + l,...,2m~l. (7) 



d?' 



(ii) There exist a + 0, s(f) G R m _„[f], and q(t) G R m _i [t] such that s(0) — q(0) — 1 and 

s 2 (t)r(t) = (af n + q(t))q(t). (8) 

Proof. We split the proof in three steps. 

Step 1: C{m,n) <=> (0). C{m,n) means that the m — n + 1 columns of the matrix given in 
Theorem|6]are linearly dependent, so there exist Sq, ■ ■ ■ , s m -„ 6 R, not all zero, such that 

so x (first column) + • ■ ■ + s„,_„ x (last column) = 0, 



which is equivalent to condition © when s{t) = Ya=o s i' G K m _„[f]. 

Step 2: (0) => (©. If s(t) e R m - n [t] verifies ©, then g(t) = X^og/?' := s(t)f(t) verifies that 
gi = for I = m + 1 , . . . , 2m - 1 . Hence, 

g(f) = 9 (r) + af/2 + 0{t 2m ), 

where = go + • • ■ + gm-\f" e R m _i[fl and a = 2g m . Therefore, 

5 2 r = s 2 / 2 = <? 2 = q 2 + at'"q + 0(t 2m ) = (q + at'")q + 0(t 2m ), 

and so s 2 r — (q + at m )q, since deg[s 2 r] < 2m - 1 and deg[(g + at m )q] < 2m - 1. Besides, a + 0, 
because deg[r] is odd and s(t) £ 0. 

Let h(t) = Z,>o hit 1 := g 2 (t) = s 2 (t)r(t) e K 2m -i M- Then 

2m 

= h 2m = ^ g;g2m-i = (gm) 2 + 2g g 2m => gogim = -a 2 /8 + 0. 

1=0 

From this property, we deduce that q(Q) = go + and s 2 (0) = q 2 (Q)/r(0) + 0. Thus, we can nor- 
malize s(f) by imposing s(Q) = 1, since condition (fTJ) only determines s(f) up to a multiplicative 
constant. This implies that q 2 (Q) = s 2 (0)r(0) = 1, so g(0) = ±1. We can assume, without loss of 
generality, that q(0) = 1 . On the contrary, we substitute q(t) and a in the identity s 2 r = (af"+q)q, 
by -q{t) and -a, respectively. 

Step 3: ((HP => (0. If there exist a + 0, s(f) e R m _„[f], and 9 (f) e R m _i[f] such that s(0) = 
q(0) = 1 and relation © holds, we set g(t) = Efeog/^ : = s (t)f(0- Then, 

g 2 = s 2 / 2 = i 2 r = (q+ at n )q = (1 + at m /q)q 2 . 
The last operation is well-defined for small values of |/|, because q(Q) + 0. Hence, 



arf" L at'" . r, m \ at 



± q ^\ + ^- = ±q \l + ^-+ o(f 2m )j = +<?+ — + 0(f 2 '"), 

so gi — for / = m + 1 , . . . , 2m - 1 . □ 

Next, we present three examples of the results that can be obtained from this formulation. 

Theorem 8. The nonsingular billiard trajectories inside the ellipsoid (0 sharing the caustics 
Qai > • • ■ i 2/t„-i are periodic with: 

• Elliptic period m — n if the roots off — Y\" = i(t — aj) are the caustic parameters; 

• Elliptic period m — n+1 if there exists Jel such that the roots off' +1 -(t—d) 2 Ylk=i ( f _/ W 
are the ellipsoidal parameters; and 

• Elliptic period m — 2n — 1 if all the roots oft 2 "' 1 — rij=i( f — a j) n't=i( f — ^k) are double. 
Proof. Suffice it to find a + 0, s(t) e R m _„[f], and e R m _i[f] such that 

s 2 (t)r(t) = (at m + q(t))q(t), s(0) = q(0) = 1, 
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for m = n, m = n + 1, and m = In - 1, respectively; see Theorem [7] We recall that r(t) = 
Til"! 1 (1 - f/ci) with {ci, . . .,c 2 „-i} = . ,«„} U {Ai, . . .,/l n _i}. 

Case m = n. If the caustic parameters are the roots of f" - I~~[" =1 - af), then 

<"-n^ia-aj) = *n£i 1 (* -/«*), 

for some factor cel. Indeed, *• = n"=i a? II^! ^ We take a = (-1)" Y\%i a] 1 , s(f) = 1, and 
q(f) = n£}(l - *Mfc). Clearly, a # 0, s(0 e KoM, 9(0 e K„-i|>L and s(0) = 9(0) = 1. Besides, 

AOKO = n "=id - t/aj) n£}(l - fMi) = (an" + 9(0)9(0. 

since Il"=i(l - t/aj) = a Il"=i(? - «;) = - a* YllZli 1 ~ fa) = + 9(0. 

Case m-n + 1. If ai, . . . ,a„ are the roots of f" +1 - (f - a") 2 n*=}( f _ fa)> men 

f +i - (f - j) 2 nd (? - ^) = * rijLiC* - «a 

for some k e E. Indeed, /c = t/ 2 n£l ^ Il"=i aj 1 - We take a = (-1)" +1 £T 2 n^j *(0 = 
(l-t/d), and 9(0 = IT/Li (!-'/<»/)• Clearly, a ^ 0, s(0 e RiW, 9(0 e E„[f], ands(0) = 9(0) = 1. 
Besides, 

n n-l 

s 2 (t)r(t) = (1 - f/af \\{\ - t/m) - t/A k ) = {af +l + q(t))q(t), 

7=1 *=1 

since (1 -r/a 1 ) 2 njEJa- W = a(?-^) 2 HEi'd -*Mt) = at" +1 -a K U n j=l (t-aj) = at n+1 +q(t). 
Case m — 2n - 1. If ? 2 "~' - I~[ 2 =i _ c d nas double roots d\,..., d n -\, then 

( i - 1 -iiH 1 (f-ci)=*rEi 1 (f-<tt 2 , 

for some k e R. Indeed, k = n^f 1 c i n? = 7 ^ 2 - We take or = - ft 2 ^ 1 c 7^ s ( f ) = IX/d ~ t/di), 
and o(f) = s 2 (0- Clearly, a + 0, s(0 e R«-iW, 9(0 6 HWzM, and s(0) = 9(0) = 1. Besides, 

Aoko = n; = 7(i - t/d,y- nSr'a - r/co = (a? 2 "- 1 + 9(0)9(0, 
since nSr^i - f/c,) = a nSr 1 ^ - c <) = at2 "~ l - aK n"=7( f - *) 2 = Q,f2 " _1 + □ 

Several questions arise about the periodic trajectories found in the previous theorem. Let 
us mention just three. Which are their caustic types, their (Cartesian) periods, and their wind- 
ing numbers? Inside what ellipsoids exist them? Are there other nonsingular periodic billiard 
trajectories with elliptic period three, four or five? 

We will give some partial answers in the next sections. 

Some technicalities become simpler after the change of variables t - l/x. Thus, we state 
another polynomial formulation of the generalized Cayley condition C(m, n). 

Proposition 9. Let R(x) — x T\ 2 =i ( x ~ where yt = 1/c,-. The generalized Cayley condi- 
tion C(m,n) holds if and only if there exist two monic polynomials S(x),R(x) e M[x] such that 
deg[5] = m — n, deg[f] = m, P(Q) + 0, and 

S 2 (x)R(x) = P(x)(P(x)-P(0)). (9) 

Furthermore, if such polynomials S (x) and P(x) exist, the following properties hold: 
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1. S(x) has no multiple roots; 

2. All the real roots ofS(x) are contained in {x € R : R(x) < 0); 

3. All the roots of S (x) are real when m < n + 3; and 

4. P(x) andP(x) — P(0) have the same number of real roots (counted with multiplicity). 

Proof. The "if and only if" follows directly from the change of variables t — l/x. Concretely, 
the relation between the objects of identities ([8} and (0 is 

R(x) = x 2 "r(l/x), S(x) = x m -"s(l/x), P(x) = a + x m q(l/x). 

Then P(Q) + if and only if a + 0, s(Q) = 1 if and only if S (x) is a monic polynomial of degree 
m - n, and q{0) = 1 if and only if P(x) is a monic polynomial of degree m. 
To prove the first two properties, suffice it to prove that gcd[S,RS'] = 1 and 

/+ := #{x e R : S(x) = < R(x)} = 0. 

If W(x) = P(x)(P(x) - P(0)) and T(x) = P(x) - P(Q)/2, we get from © that 

W(x) = S 2 (x)R(x) = T 2 (x) - P 2 (0)/4, 
W'(x) = S (x)(S (x)R'(x) + 2R(x)S '(x)) = 2T(x)P'(x). 

We consider the factorization W'(x) = 2mW-(x)Wo(x)W + (x)W*(x), where if z e C is a root 
of multiplicity y3 of W'(x) such that W(z) < 0, W(z) = 0, W(z) > 0, or W(z) t K, then (x - zf 
is included in the monic factor W-(x), Wq(x), W+(x), or W t (x), respectively. Next, we find some 
lower bounds of the degrees of these factors. 

First, T is divisor of W_, because W takes the negative value -P 2 (0)/4 at each root of T. 
Hence, deg[W_] > deg[T] - m. Second, S gcd[S,RS'] is a divisor of Wo, because W vanishes at 
each root of S . Thus, deg[Wo] > m - n + Iq, where Iq denotes the degree of gcd[5, RS ']. Third, 

C W(a) = W(b) = 1 

deg[W + ] >#<(a,b)<zR: R(x) > for all x e (a, b) i=n-l+l + . 
[ S(x) *0forallxe (a,b) j 

To understand the above inequality, we realize that if (a, b) is an open bounded interval that 
satisfies the above three properties, then W(x) = S 2 (x)R(x) > for all x € (a,b), and W'(x) 
vanishes at some point c e (a, b), by Rolle's Theorem. Therefore, deg[W+] is at least the number 
of such intervals. We combine these three lower bounds: 

2m - 1 = deg[W] > deg[W_] + deg[W ] + deg[W + ] > 2m - 1 + l + l + . 

This implies that / = 1+ = 0. Indeed, W_ = T, W = S, W, = 1, and gcd[S,RS'] = 1. 

Next, we prove the property about the number of roots of P(x) and P(x) - P{0). Let z be 
a root of the derivative P'. Since W = 2TP' and W- = T, we deduce that W(z) cannot be a 
negative number. This implies that if P(z) is a real value between and P{0), then P'(z) + 0, 
since W{z) = P(z)(P(z) - P(0)) < 0. In particular, we deduce that the number of real roots 
(counted with multiplicity) of the polynomial P(x) - rj does not change when the constant 77 e R 
moves from to P(0). 

Finally, we prove that S (x) has only real roots when m < n + 3. Let us suppose that z t R is a 
root of S (x). Then z is also a root of S (x), so (x - z)(x -z)\S (x). Using the identity S 2 (x)R(x) = 
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P(x)(P(x) - -P(O)), we get that (x - z) 2 (x - z) 2 is either a divisor of P(x) or a divisor of P(x) - P(Q), 
since P(x) and P(x) - P(Q) have no common factors. But P(x) and P{x) - P(Q) have the same 
number of real roots, so there exists another w t K U {z, z] such that (x — w) 2 (x — w) 2 is a divisor 
of P(x)(P(x) - P(0)). This implies that S(x) has, at least, four different complex roots, and so 



There are some theoretical arguments against the existence of non-real roots of polynomial 
S (x), although we have not been able to prove it. 

Conjecture 2. Let R(x) = x Tl^^i* ~ 7i) with < yin-i < ■ ■ • < J\- If relation (0 holds for 
some polynomials S (x), P(x) e R[x] such that P(0) + 0, then S (x) has only real roots. 

5. Generalized Cayley conditions in the minimal case 

Let us consider the case of minimal elliptic periods; that is, m-n. 

We begin with a technical lemma to describe how the roots of the polynomials of the form 
P(x)(P(x) - P(0)) with P(Q) + are ordered in the real line, assuming that all these roots are 
positive — but the trivial one — , and have multiplicity at most two. 

Lemma 10. Let P(x) e be a monic polynomial of degree m such that P(0) + and all the 
roots of P(x)(P(x) — P(0)) are positive — but a simple root at x — — , and have multiplicity at 
most two. Let a m < ■ • • < a>\ be the positive roots of P(x). Let fi m -\ < ■ ■ ■ < /3\ be the positive 
roots ofP(x) - P(0). 

If m is odd, then f52i-\,Pii £ (0:21, 021-1) for all I = 1, . . . , (m — l)/2; so 

< a m < ff,„_i < j8„,_i < p m -2 < a m - 2 < a m - 3 < ■ ■ ■ < a 3 < a 2 < P2 < Pi < a\. 

If m is even, then f5\ > a\ and j$2ufi2h\ ^ (o<2;+i, 0:21) for all I — 1, . . . , (m — 2)/2; so 

< a m < < p m -\ < p m -2 < a m - 2 < a,„-3 < ■ ■ ■ < < /3 2 < a 2 < a x < /3i . 

Proof. Let 77 e K. Using that the only critical points of P(x) are non-degenerate local maxima or 
non-degenerate local minima, we deduce that the polynomial P(x) - r\ has m real roots (counted 
with multiplicity) if and only if r\ < 77 < rj, where 



Therefore, i] < min(0, P(0)) and i] > max(0, P(0)). 

We begin with the case m odd, so P(0) = Y\%\ «j < 0, 77 < P(0), and ~rj > 0. The 

roots of P(x) and P(x) - P(Q) can be viewed as the abscissae of the intersections of the graph 
{y = P(x)} with the horizontal line {y = 0} and {y - P(0)}, respectively. Double roots correspond 
to tangential intersections. We know that P(x) > rj > at the local maxima, and P(x) <rj< P(Q) 
at the local minima. This means that the intersections of the graph {y = P(x)} with the lines 
{y = 0} and {y = P{0)} have the following pattern from left to right. First, the graph crosses 
{y = P(Q)} at the abscissa x = 0; second, it intersects {y = 0} at two abscissae a m and ar m _i, 
which may coincide giving rise to a double root of P(x); third, it intersects {y = P(Q)} at two 
abscissae /3 m _i and /3,„_2, which may coincide giving rise to a double root of P(x) - P(0); fourth, 
it intersects {y = 0} at two abscissae o , m _2 and a,„_3, which may coincide giving rise to a double 
root of P(x); and so on. The last intersection correspond to the abscissa x = a\ . 

The proof for m even is similar. We skip the details. □ 



m-n — deg[5] > 4. 



□ 
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We emphasize that ellipsoidal parameters < a\ < ■ ■ ■ < a„ and nonsingular caustic pa- 
rameters A\ < ■ ■ ■ < /l„_i verify restrictions (0; then the parameters < c\ < ■ ■ ■ < C2n-\ are 
defined in (0; next the inverse quantities < jin-\ < ■ ■ ■ <j\ are given by y, = 1/c,; and finally, 
R( x ) = x Uil l (x - yd. We will make use of these notations, orderings, and conventions along 
the paper without any explicit mention. 

Corollary 11. Let {1, . . ,,2n — 1} = /„ U K„ be the decomposition definedby 

/i={l}, J 2 = {2,3}, y„ = 7„- 2 U{2«-2,2n-l}, K n = J n . x . 
IfP(x) e R[x] is a monic polynomial of degree n such that P(0) + and 

R(x)=P(x)(P(x)-P(0)), 
then P(x) = UjeJ„( x ~ 7j) = P(P) + x UkeK„(x - y k ). 

Proof. There exists a decomposition {1, . . . ,2n — 1} = J' UK' such that #/' = n, #K' = n — 1, 
and P(x) = Y\jer( x ~ 7j) = ^(0) + x I\keK'( x ~ 7k)- The polynomial P(x) verifies the hypotheses 
stated in Lemma [TOl so the roots {a\, . . . , a m ) = {jj : j e J'} and {/3i, . . . ,p m -\) = {7k '■ k s K') 
obey the ordering described in that lemma. Therefore, J' = J„ and K' — K„. □ 

We rewrite now the generalized Cayley condition C(n, n) using the previous results. For 
brevity, we omit the dependence of the decomposition {1 , . . . , 2n — 1} = J UK on the index n. We 
note that #J = n and #K -n-l. The symbol e/("a set of parameters") denotes the elementary 
symmetric polynomial of degree I in those parameters. 

Proposition 12. C(n, n) is equivalent to any of the following four conditions: 

(i) IfP(x) = W je] {x - jj), then P(x) - P(Q) = x Y\ keK (x - y k ). 

(ii) Ujej(7j ~ 7k) = Ujej7j, far all k e K. 

(iii) &i({jj}jej) = zi{{7k\keK), for all l=l,...,n-l. 

(iv) Zjcj 7 l j = ZkcK 7 l k , for all 1=1,..., n-l. 

Proof. We split the proof in four steps. 

Step 1: C(n, n) <=> (0. Let us assume that C(«, n) holds. Then there exist a monic polynomial 
P(x) e R[x] of degree n such that P(0) + and 

R(x) = P(x)(P(x) - P(0)). 

Thus, condition (H) follows from Corollary QT| 

Reciprocally, if condition © holds, P(x)(P(x) - P(0)) = x Yii^ix ~ so C(n, n) holds. 
Step 2: (0 o (O. If P(x) = \\ je] (x - yj) and Q(x) = x Ukexix - 7k), then 

Q(x) = P(x) - P(0) o P(y k ) = P(0), Vfe e K o UjeAjj ~ 7k) = Ujej Yj, e K. 

Step 3: (0 o If P(x) = Y\jeA x ~ Yj) and Q( x ) = x UkeK ( x ~ 7k), then 

n— 1 n— 1 

Q(x) =x" + Yj- l ) l ^7k)keK) xn - 1 , P(x) = x" + J](-l) / e / ({ r; } iey )x"-' + P(0). 
i=\ i=i 

Step 4: ( 17771 ) o (jTv}. It follows from the Newton's identities connecting the elementary sym- 
metric polynomials and the power sum symmetric polynomials; see j26ll . □ 
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Example 1. The quantities 73 = 1, 72 = 2, and j\ — 3 verify C(2, 2), because 1+2 = 3. This 
means that the billiard trajectories 

• Inside the ellipse Q : x 2 + 2y 2 = 1 with caustic parameter A = 1/3; or 

• Inside the ellipse Q : x 2 + 3y 2 = 1 with caustic parameter A = 1/2; 

are periodic with elliptic periodic m-2. Their caustic types are g — and g = 1, respectively. 

Example 2. The quantities 75 = 1, 74 = 2, 73 = 4, 72 = 5, and 71 = 6 verify C(3, 3), because 
1+2 + 6 = 4 + 5 and l 2 + 2 2 + 6 2 = 4 2 + 5 2 . Hence, the billiard trajectories 

• Inside the ellipsoid Q : x 2 + 2y 2 + 5z 2 = 1 with caustic parameters A\ = i and A.% = |; or 

• Inside the ellipsoid : x 2 + 2y 2 + 6z 2 = 1 with caustic parameters ^1 = 5 and A.2 = j, or 

• Inside the ellipsoid Q : x 2 + 4y 2 + 5z 2 = 1 with caustic parameters Ai - I and A2 = 5; or 

• Inside the ellipsoid 2 : x 2 + 4j 2 + 6z 2 = 1 with caustic parameters ^1 = 5 and A^-\; 

are periodic with elliptic periodic m - 3. Their caustic types are g = (0, 1), ^ = (1, 1), §■ = (0, 2), 
and §• = (1,2), respectively. We recall that these caustic types were denoted EH1, H1H1, EH2, 
and H1H2 in the introduction. 

Let us compare the system of homogeneous symmetric polynomial equations (0, which was 
obtained directly from the matrix formulation, with the system of homogeneous non-symmetric 
polynomial equations Yjjej Tj - Tiketc Yk>l < I <n — l, obtained in the previous proposition. We 
are dealing with the case m = n, so the /-th equation of the former system has degree n+l, whereas 
the Z-th equation of the new system has degree I. Besides, the new system has a remarkably simple 
closed expression. This shows that the polynomial formulation simplifies the problem. 

The beauty of the conditions regarding the elementary symmetric polynomials and the power 
sum symmetric polynomials given in Proposition[l2]has been the motivation for the introduction 
of the inverse quantities 7, = 1/c,. Nevertheless, we find useful to state the following result in 
terms of the ellipsoidal parameters aj, in order to answer some questions about the nonsingular 
periodic billiard trajectories found in the first item of Theorem[8] 

Theorem 13. There exist nonsingular periodic billiard trajectories inside the ellipsoid (0 with 
elliptic period m — n and caustic type 

_{ (1, 1, 3,3, . . .,n — 2,n — 2) for n odd 
g ~ \ (0,2,2,4,4,. ..,n-2,n-2) forneven { ' 

if and only if all the roots of f — n}=i( f — a j) are rea ^ an d simple. These periodic billiard 
trajectories have the roots oft" — ]^[" =1 (? — aj) as caustic parameters, period mo — 2n, and even 
winding numbers mo, . . ., m n -\. Indeed, 

mj — 2m j — 2(n — j), j = 0, . . . ,n — 1, (11) 

provided Conjecture\l\on the strict decreasing ordering of winding numbers holds. 

Proof. Let us assume that there exist nonsingular periodic billiard trajectories with elliptic period 
n and caustic type (TTOb . By definition of caustic type, the caustic parameters A\ < ■ ■ ■ < A„-\ of 
such trajectories verify that 
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• If n is odd, then A,2i-\, A21 e (0.21-1,021), for / = 1, . . . , (n — l)/2; 



• If n is even, then A\ e (0, a\), and A21, /I2/+1 e (021, 021+1), for Z = 1, ...,(«- 2)/2. 

Hence, we can split the set {y,- = 1/c; : i = 1, . . . , 2n - 1} as the disjoint union of the sets 

{yy : ; e J) = [l/a u l/a„), \y k : k e K) = {1/*,. . . , 

where {1, . . . , 2n - 1} = J U K is the decomposition described in Corollary [TT] Thus, we know 
from condition|Iv]of PropositionfT2lthat 

frfl_ ILfrl, jfc=i,. ..,„-i. 
1=}U (=}«/ 

This identity can be written as A' l k = rifciWfc _ a j)< f° r all & = 1, 1, which implies that 
the caustic parameters Aj, . . . , A n -i are the roots of f — ITj=i( f _ a j)- 

Reciprocally, let us assume that the roots of q(t) = t" - ]~7" =1 (f - aj) are real and simple. Let 
A\ < ■ ■ ■ < An-i be these roots. None of them is zero, since g(0) + 0. Besides, Ak is a root of q(t) 
if and only if := I /At + is a root of 

eW;= (-')-v g ( 1M = pw _ p(0)| 

^/i ■ ■ • 

where = EI^iC^ _ a j) w i tri ff i = ^l a j- Therefore, the roots aj = 1/aj and fi^ = 1/Ak are or- 
dered as stated in LemmafTUl The consequences are two-fold. On the one hand, Ak € (a ft , a ft+ i), 
where ^ = (^1, . . . ,g„-\) is the caustic type given in ( TTOb . On the other hand, there exist non- 
singular billiard trajectories inside the ellipsoid Q sharing the caustics Qx x , . . . , Qa„_ x , because 
the existence conditions (f2j) hold. Thus, the trajectories sharing the caustics Q^, . . . , Qa„^ are 
periodic with elliptic period n and caustic type g, since the generalized Cayley condition C(n, n) 
holds; see PropositionfL?! 

Next, we prove the claims on the (Cartesian) period and the winding numbers. The caus- 
tic parameters are located in the intervals delimited by the ellipsoidal parameters given at the 
beginning of the proof, which implies that 

\c2j,c2j+\) n . ..,«„} t 0, j = 0,...,n- 1, 

where cq := < c\ < ■ ■ ■ < C2 n -\ are defined in (f3]). Thus, all winding numbers are even — see 
Theorem|5] — and so, by definition of elliptic period, mo = 2m = 2n. 

Finally, let us assume that winding numbers are ordered as stated in Conjecture Q] so 2 < 
m„_i < • • • < mo — 2n with too, . . . , to„_i even. Then mj — 2rhj — 2{n — j). □ 

Remark 1. If all the roots of f - Y\" =x (t - aj) are real, but some of them are double, then we get 
singular periodic billiard trajectories. In that case, there are only two possible scenarios. Either 
n is odd and A%-\ = A% for some I = 1, . . . , (n - l)/2; or n is even and A21 = /I2/+1 for some 
/ = 1, . . . , (n - 2)/2. In all these cases, the singular periodic trajectories are formed by segments 
contained in some nonsingular ruled confocal quadrics. 

Remark 2. All the periodic billiard trajectories mentioned in Theorem[13]have caustic type ( TTOb . 
One may establish similar theorems for other caustic types. For instance, the versions EH1, 
H1H1, EH2, and H1H2 of Theorem[T3lin the spatial case will be listed in Table|2] 
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6. Cayley conditions in the general case 



Once we have understood the minimal case m — n, we tackle the general case to > n. 

Let us explain the fundamental question by means of an example. In Section [3] we saw that 
condition C(3, 2) becomes a single homogeneous symmetric polynomial equation of degree eight 
in the variables 71,72,73, with only two feasible solutions; namely, the ones given in (|6). Thus, 
it is natural to ask whether can we rewrite C(3,2) as a set of two simpler conditions such that 
each one of them gives rise to one of the solutions given in ©. 

By the way, we raise this question for any m > n. Can we rewrite C(m, ri) as a set of "simpler" 
conditions such that each one of them gives rise to just "one" solution of C(m, ri)l We answer this 
question in the affirmative. Indeed, we parametrize these "simpler" conditions by the elements 
of the set 

T(m, ri) — J(ti , . . . , t„) eZ" : t\ + • • ■ + t„ — m - n, ti,...,t„ > o}. 

The cardinal of T{m, ri) is the number of monomials of degree m - n in n variables. Thus, 
#T (to, ri) = ulj), which gives a precise estimate of the complexity of the Cayley condition 
C(m, ri) when to grows. We will refer to the elements of T(m, ri) as signatures. We set 72,, — in 
order to simplify some notations. 

Definition 5. Given any signature t = (ti , . . . , t„) e T(m, ri), we say that condition C(m, n; t) 
holds if and only if there exist two monic polynomials S (x), P(x) e M.[x] such that degfS ] — m-n, 
deg[P] = to, P(0) + 0, S 2 (x)R(x) = P(x)(P(x) - P(0)), and S(x) has to - n simple real roots 
<5 m -« < • • ■ < 61 such that 

#({<5i,...,<5 m _J n (72r,72r-i)) =i" r , r=l,...,n. (12) 

Corollary 14. If there exists t e Tim, ri) such that C(m, n; t) holds, then C(m, ri) also holds. The 
reciprocal implication is true for to < n + 3 (or provided Conjecture\2\holds). 

Proof. The first implication is obvious. For the reciprocal implication, we simply recall that S (x) 
has only real roots when to < n + 3 and all its real roots are contained in \x e R : R(x) < 0} = 
U"=i(72r, 72r-i); see Proposition [9] □ 

Definition 6. Given any signature t € T(m, ri), let {1 , . . . , 2n — 1} = J T U K T and { 1, . . . , m — n] = 
V T U W T be the decompositions determined as follows. If 6 m -„ < ■ ■ ■ < S\ is any ordered sequence 
verifying ( fT2l . then the elements of the multisets 

{ai,..., a m } = [yj: j€j T }\J{6 v ,6 v : ve V T }, 
{jSi, • ■ ■ ,j8 m _i } = {jk'-ke K T ] U {6 W , S w : w e W T ], 

are ordered as in Lemma [TUl 

Multisets are a generalization of sets in which members are allowed to appear more than 



once; see H27H . In our case, the numbers Si,..., S,„-„ appear twice. 

These decompositions are well-defined. That is, they only depend on the signature t, since 
any ordered sequence S,„-„ < ■ ■ ■ < Si verifying (TT2l gives rise to the same decomposition. The 
decomposition { 1 , . . . , 2n - 1 } = J„ U K„ given in Corollary QT]correspond t° the trivial signature 
: = (0 0)Ef(n,4 

Next, we generalize Corollary QT| and Proposition[T2]to the case to > n. 
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Corollary 15. Let 

Sm-n < ' ' ' ^ &i b€ an ordered sequence verifying ( I72D for some signature 
r € T"(ra, n). IfP(x) is a monic polynomial of degree m such that P(0) and 

U::"(x - 5 U ) 2 ■ R(x) = P(x)(P(x) - P(0)), 

then P(x) = UjeJ T ( x ~ Yj) YlvsvM ~ 5 = P (°) + x UkeK T ( x - Jk) YlwewM ~ 5 »') 2 - 

Proof. There exist two decompositions {1, . . . ,2n - 1} = J' UK' and {I, , . . ,m — n] = V UW' 
such that P(x) = UjeAx - Yj) IlvevC* - S v f = P(0) + x rW (* - Yk) Uwewix - 6 W ) 2 . The 
polynomial P(x) verifies the hypotheses stated in Lemma[10l so the roots 

{ai,..., a m ] = {yj : j e J'} U {6 V , 6 V : v € V'} 
{flu ■ ■ ■ ,Pm-i\ = iYk : k e K'} U {6 W , S w : w e W'} 

obey the ordering described in that lemma. Hence, J' = J T , K' = K T , V = V T , and W = W T . □ 

Proposition 16. Condition C(m, n; t) holds if and only if there exists a sequence 5 m - n < ■ ■ ■ < 5\ 
verifying ( 1721 ) such that the following three equivalent properties hold: 

(i) IfP(x) = UjejM ~ Yj) UvevM ~ S *) 2 > then P ( X ) ~ P (°) = X UkeK T ( x ~ Yk) W w ewM ~ 6 -) 2 - 

(ii) &i{{jj}jej T U {S v ,6 v } v€Vt J = ei({yk)keK T U [6 W , 5 w } weWt J, for all I = l,...,m- 1. 

(iii) lLjej r y'j + 2 EvsV T K = ZkcK T Y' k +2 lweW T for all I = I, . . . ,m - I. 

Proof. We simply repeat the steps of the proof of Proposition[12] but using Corollary [T5linstead 
of Corollary E] □ 

Example 3. The quantities 73 = 1, ji = 4, and j\ — 9 verify condition C(3, 2; t) with r = (1,0), 
because 1 + 4 + 9 = 2 ■ 7, l 2 + 4 2 + 9 2 = 2 ■ 7 2 , and 7 e (4, 9). Hence, the billiard trajectories 

• Inside the ellipse Q : x 2 + Ay 1 = 1 with caustic parameter A = 1/9; or 

• Inside the ellipse Q : x 2 + 9y 2 = 1 with caustic parameter A = 1/4; 

are periodic with elliptic periodic m — 3. Their caustic types are g = and g = 1, respectively. 

All conditions C(m, n; r), r e 7~(m, «), give rise to nonsingular periodic billiard trajectories 
with elliptic period m, so we wondered which is the dynamical meaning of the signature r. 
We believe that there exists a one-to-one correspondence between the elliptic winding numbers 
mo, . . . , m„_i — see Definition[3] — and the signature r = (ti, . . . , t„). 

Conjecture 3. Set rh n — 0. Then thj — mj + \ + Tj+\ + 1 for all j = 0, . . . , n — 1. 

This conjecture follows from the interpretation of C(m, n; r) as a singular limit of C(m, m) 
when m - n couples of simple roots collide, so they become double roots. Unfortunately, we 
have not been able to transform this argument into a rigorous proof, although all our analytical 
and numerical computations agree with the conjecture. 

To end this section, we stress that if conjectures|2]and[3]hold, then the elliptic winding num- 
bers mo, . . . , m„_i of any nonsingular periodic billiard trajectory verify the above-mentioned rela- 
tions for some signature r = (ti, . . . , t„) with non-negative entries, so the sequence mo, ... , m n -\ 
strictly decreases, and ConjectureQ] holds. 
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7. The planar case 



We adapt the previous setting of billiards inside ellipsoids of W to the planar case n — 2. To 
follow traditional conventions in the literature, we write the ellipse as 

Q = l( x ,y) e R 2 : — + — = l), a > b > 0. (13) 
[ a b ) 

Any nonsingular billiard trajectory inside Q is tangent to one confocal caustic 

e, = {(^)e R -^ + i = i}, 

where AeA = EuH, with E = (0, b) and H = (b, a). 

The names of the connected components of A come from the fact that Qx is a confocal ellipse 
for A e E and a confocal hyperbola for A e H. The singular cases A - b and A - a correspond 
to the x-axis and y-axis, respectively. We say that the caustic type of a billiard trajectory is E or 
H, when its caustic is an ellipse or a hyperbola (compare with Definition [l]). We also distinguish 
between E-caustics and H-caustics. 

We recall some concepts related to periodic trajectories of billiards inside ellipses. These 
results can be found, for instance, in 1 91 12411 . To begin with, we introduce the function p : A — > W. 
given by the quotient of elliptic integrals 



,-min(M) i t 




p(A) =p(A;b,a) := rT ^ . (14) 



Jmax(M) ^(A-t)(b-t)(a-l) 

It is called the rotation number and characterizes the caustic parameters that give rise to periodic 
trajectories. To be precise, the billiard trajectories with caustic Qx are periodic if and only if 

p(A) — mi/2mo e Q 

for some integers 2 < m\ < niQ, which are the winding numbers. On the one hand, mo is the 
period, On the other hand, m\ is twice the number of turns around the ellipse Qx for E-caustics, 
and the number of crossings of the y-axis for H-caustics. Thus, mi is always even. Besides, all 
periodic trajectories with H-caustics have even period. (Compare with Theorem[5]) 

Proposition 17. The winding numbers 2 < m\ < mo, rotation number p — m\/2tno, signature 
t — (ti , T2) € T(m, 2), caustic type (E or H), and caustic parameter A of all nonsingular periodic 
billiard trajectories inside the ellipse ( liiD with elliptic period m e {2,3} are listed in Table\l\ 
The ellipses where such trajectories take place are also listed. 



Proof. We split the proof in four steps. 

Step 1: To find the solutions ofC(m, 2) in terms of the inverse quantities y,-. First, we saw in 
Proposition [T2l that C(2, 2) holds if and only if j\ =72+73- 

Next, we focus on the case m — 3. We note that C(3, 2) holds if and only if C(3, 2; t) holds 
for some t = (ti, T2) e Z 2 such that t\ + T2 = 1 and n, T2 > 0; see Corollary fT4l 



18 



m 


mo 


mi 


P 


T 


Type 


Ellipses 


2 


4 


2 


1/4 


(0,0) 


E 


any 


2 


4 


2 


1/4 


(0,0) 


H 


2b < a 


3 


6 


2 


1/6 


(1,0) 


E 


any 


3 


6 


2 


1/6 


(1,0) 


H 


4b < a 


3 


3 


2 


1/3 


(0,1) 


E 


any 


3 


6 


4 


1/3 


(0,1) 


H 


4b < 3a 



Caustic parameter 

ab 

«+/> 
ab 

ab 



a+b+2 Voi 
ab 

a+b-2 slab 
3ah 



a+b+2 ia 2 -ab+b 2 
ab 

2^a 2 -ab+b-a 



Table 1: Algebraic formulas for the caustic parameter corresponding to nonsingular periodic billiard trajectories with 
elliptic period me {2, 31 in the planar case. 



Let us begin with the signature t = (1,0). After a straightforward check, we get that the 
decompositions presented in Definition [6] are J T = {1,2,3}, K T — V T — 0, and W T = {1}. Thus, 
C(3, 2; t) holds if and only if there exists some 6i e (72, 7i) such that 

P(x) = (x - yi)(x - y 2 )(x - y 3 ) = P(0) + x{x - 6 x f, 

or, equivalently, if and only if the discriminant of the polynomial 

P(x) — P(0) 

Q(x) = = x 2 - ei (yi,y 2 , 73)* + e2(7i>72,73) 

x 

is equal to zero. The discriminant of Q(x) is 

A = y\ + y\ + y\ - 2y\y 2 ~ 2ym ~ lyiyi- 

We already saw in Section[3]that the only feasible solution of A = is -Jyi = -^72 + y[Ti- 

When t = (0, 1), the decompositions are J T = {1}, K T = {2, 3), V T = {1}, and W T = 0. Thus, 
C(3, 2; t) holds if and only if there exists some 61 e (0, 73) such that 

71+2^1=72+73, 7^ + 2(5; = 72 + 73, 
or, equivalently, if and only if 

37? - 2 (72 + 73)71 - (72 - 73) 2 = (72 +73 -7i ) 2 -2(7 2 + 73 _ 7i) 

= (25i) 2 - 4<^ = 0. 

And we already saw in Section|3]that the only feasible solution of the above equation is 

37i = 72 + 73 + 2 yjy\ +y\- y 2 y^. 

Step 2: To express the above solutions in terms of a, b, and A. If the caustic type is E, then 
A e (0,b), yi = 1 1 A, y 2 = l/b, and 73 = I /a. Thus, 

7i = 72 + 73 

VtT = V72 + V73 

371 =72+73+2 yjyl + yj- 7273 
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^ n — ab 



a+b ' 
ab 

a+b+2 yfab 
3ab 



a+b+2 ^a 2 -ab+b 2 ' 



Figure 1: Periodic trajectories corresponding to the caustic parameters listed in Table[T] The ellipse for A = ^°^^ f — is 
flatter, because it must satisfy condition 4b < a. 



If the caustic type is H, then A e (b, a), y\ — l/b, j2 — 1 I A, and 73 = 1 I a. Thus, 



a-b' 
ab 



3yi = 72 + 73 + 2 >/r| +7^-7273 » ^ = 



2 \'a 2 -ab+b-a ' 



Step 3: To determine the ellipses where such periodic billiard trajectories take place. We ask 
whether the caustic parameters found above belong to the interval (0, b) for E-caustics, and to 
the interval (b, a) for H-caustics. The caustic type E does not give any restriction, because 

< b < a and A e f-^-, — &-=,, 1===} => A e (O.fc). 

^a+fc' a+b+2yfab' a+b+2y/a 2 -ab+b 2 j v ' 

On the contrary, the caustic type H gives rise to some restrictions. Namely, 

b < ^ < a <=> 2b < a, 
b < — < a o 4fr < a, 

a+b-2 yab 

b < —f= < a o 4b < 3a. 

2 ya 2 -ab+b-a 

Step 4: To find the winding numbers and the rotation number. The winding numbers 2 < 
m\ < mo and the rotation number p{A) = oti/2oto are obtained from geometric arguments. To be 
precise, we draw in Figure[T]a billiard trajectory tangent to Q.i for each of the caustic parameters 
listed in Table[T] Then we recall that mo is the period and m\ is twice the number of turns around 
the ellipse Q,\ for E-caustics, and the number of crossings of the y-axis for H-caustics. □ 
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Let p* €{1/3,1 /4, 1 /6). We have seen above that inside any ellipse dT3l > there exists a unique 
E-caustic whose tangent billiard trajectories have rotation number p». Besides, if A E (a, b; p*) and 
Au(a, b; p*) denote the caustic parameters associated to the E-caustic and H-caustic with rotation 
number p», we see that 

A E {b,a\p*) = A E (a,b;p*), b = A E (a, A H (a, fc;p»);p„). (15) 

These properties can be generalized. On the one hand, the rotation number (TT41 is an in- 
creasing function in the interval (0, b) such that p(0) = and p(b) = 1/2; see ll24ll . This means 
that given any p„ € (0, 1 /2), there exists a unique A* e (0, b) such that p(/U) = p«, and so, there 
exists a unique E-caustic whose tangent billiard trajectories have rotation number p». On the 
other hand, relations ( fTBI ) can be obtained by using that the rotation number (fT~4-b is symmetric 
in the three parameters a, b, and A. Consequently, one can find the formula for An (respectively, 
A E ) from the formula for A% (respectively, Ah) by using the second relation. 
It is interesting to realize that the results in Table[T]agree with Conjecture[3] 
In the planar case n = 2, the caustic type ( fTUb is g — or, equivalently, E. Hence, the planar 
version of TheoremfOlis shown in the first row of TableQ] because A = A E (a, b; 1/4) = ab/(a+b) 
is the root of f 2 - (f - a)(t - b). This naive observation was the germ of this paper. 

8. The spatial case 

In order to study the spatial case n — 3, we consider the triaxial ellipsoid 

Q = Ux,y,z) 6 K 3 :— + — + — = l), a>b>c>0. (16) 
[ a b c ) 

Any nonsingular billiard trajectory inside Q is tangent to two distinct nonsingular caustics 
Qaj and Q^, with Ai < A2, of the confocal family 

ai = f(*,y,z)eR 3 :^ + -^- + ^- = l). (17) 
( a — Ab — Ac — A J 

The caustic Q \ is an ellipsoid for A e (0, c), a hyperboloid of one sheet when A e (c, b), and a 
hyperboloid of two sheets if A e (b, a). Not all combinations of nonsingular caustics can take 
place, but only the four caustic types EH1, H1H1, EH2, and H1H2. 

Proposition 18. The caustic type and caustic parameters of all nonsingular periodic billiard 
trajectories inside the triaxial ellipsoid A16\l with elliptic period m — 3 are listed in Table\2\ The 
ellipsoids where such trajectories take place are also listed. 

Proof. If a, b, and c are the ellipsoidal parameters, and A\ and A2 are the caustic parameters, we 
set {c\,C2, C3, C4, C5} = {a, b, c, A\, A2}, where < c\ < C2 < C3 < C4 < C5. We also set y t = 1/c/. 
Let {1,2,3,4,5} = J U K, with J = {1,4,5} and K = {2,3}, be the decomposition defined in 
Corollary QT| when n — 3. From Proposition[12]we know that 

C(3, 3) o y 2 +73 = 7i + 74 + 75 and y\ + y\ = y\ + y\ + y] 
<=> (71 - 7/0(74 - 7/0(75 - Jk) = 717475, for k = 2, 3 
o c\ = (q - d)(c 4 - Q)(C5 - a), for = 2,3. 
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Type Ellipsoids 



Caustic parameters 



EH1 
H1H1 

EH2 
H1H2 



ab 



c < 



c < 



a + b + yfab 
ab 

a + b + 2 ^fab 
a -2b 

c < — a 

2a — 3b 
2b < a 

a — 2b 
c < -ab 



b > 



(a - b) 2 
ac 



a + c — yfac 



c 3 — (c — A\)(b — c)(a — c) 
l/A 2 + l/c= l/a+l/b+ l/A l 

Roots of f 3 - (? - a)(t - b)(t - c) 



Roots of (a - b)(a - c)t + (be - a(b + c))at + a be 



b 3 = (b- c)(A 2 - b)(a - b) 
l//li + l/i= l/a + l/c+l/A 2 



Table 2: Algebraic formulas for the caustic parameters corresponding to the nonsingular periodic trajectories with elliptic 
period m = 3 in the spatial case. 



In the rest of the proof, we study each caustic type separately. 
Caustic type EH1. In this case < A\ < c < A 2 < b < a, so 

c\ — A\, c 2 — c, C3 = A 2 , C4 — b, C5 = a. 

Thus the formula for A\ follows from relation = (c 2 - C\)(c4 - c 2 )(c$ - c 2 ), whereas the formula 
for A 2 follows from relation 72+73 = J\ + J a + 75 ■ Next, we look for ellipsoidal parameters such 
that the caustic parameters computed using these two formulas are placed in the right intervals: 
A\ e (0, c) and A 2 e (c, b). 

To begin with, we note that A\ < c, since (c - Ai){b - c)(a - c) — c 3 > 0. Besides, 

ab — (a + b)c ab 
Ai = c > o c < 



(b - c)(a - c) a + b 

On the other hand, if A\ € (0, c), then 

1/A 2 = + + - 1/c) > 1/a+l/b > \/b, 

so A 2 < b. Finally, 

1 c 111112 11 
A 2 > c <=> - + — — = — = — + < - 

c ab — (a + b)c A\ A 2 c a b c a b 

ab 

o c < 



a + b + yfab 



Therefore, A\ e (0, c) and A 2 e (c, b) if and only if c < ab/(a + b + yfab). 

Caustic type H1H1. If n — 3, then the caustic type ( flOb is g — (1, 1) or, equivalently, H1H1. 
Hence, the study for the caustic type H1H1 was already carried out in TheoremQj] Suffice it to 
note that the polynomial 

f 3 - (f - a)(t - b)(t -c) - (a + b + c)t 2 - {ab + ac + bc)t + abc 

has two real simple roots if and only if its discriminant 

A = (ab + ac + be) 2 - 4abc(a + b + c) - (a - b) 2 c 2 - 2ab(a + b)c + a z b 2 
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is positive. This discriminant is a second-degree polynomial in c whose roots are 

ab(a + b) + 2ab y/ab ab 

(a-b) 2 a + b + 2yfab' 

We note that < c_ < b < c+. Thus, using that < c < b < a, we get A > <=> c < c_. 
Caustic type EH2. In this case < A\ < c < b < A2 < a, so 

C\ = X\, c 2 — c, C3 — b, C4 — A2, C5 = a. 

Using relations y' 2 + 7^ = 7' +y l 4 +y l 5 , with / = 1, 2, we know that 

11111 

Si := —r + — -r = -7 H : :, I — 1,2. 

A[ A' 2 c> b l a 1 

Hence, 1 jA\ and 1 jAi are the roots of the polynomial 

2 a?-J2 2 fcc - q(fc + c) (q - fc)(q - c) 
(x- l/ii)(x- 1/A 2 ) = x -six+ — - — =x + ■ x+ . 

2 abc a^bc 

Thus, using the change of variables t = l/x, we get that Ai and /I2 are the roots of 

Q(t) ={a- b)(a - c)t 2 + (be - a(b + c))at + a 2 bc. 

We look for ellipsoidal parameters such that Q(t) has a root in (0, c) and a root in (b, a). The 
root in (0, c) always exists, since Q(0) = a 2 bc > and Q(c) = -c 3 (a - b) < 0. Besides, 
Q(b) = -b 3 (a - c) < and lim ( ^ +0O Q(t) = +00, so g(f) has a root in (fo, a) if and only if 



Q(a) = a 2 [a 2 - 2a(b + c) + 3bc) > 0, 



or, equivalently, if and only if c < (a-2b)a/(2a-3b) and 2b < a. We have used that < c < b < a 
in the last equivalence. 

Caustic type H1H2. In this case < c < A\ < b < A2 < a, so 

C\ — c, C2 — A\, C3 = b, C4 = A2, C5 = a. 

Thus the formula for A2 follows from relation c 2 = (c 2 - c\)(ca — c 2 )(c5 - c 2 ), whereas the formula 
for A\ follows from relation 72+73 =71+74 + 75- Next, we look for conditions on the ellipsoidal 
parameters such that the caustic parameters computed from the previous formulas are placed in 
the right intervals: A\ e (c, b) and A 2 e (b, a). 

To begin with, we note that A 2 > b, because (b - c)(A 2 - b)(a - b) = b 3 > 0. Besides, 

(a + c)b - ac , a - 2b , 

— -b = A 2 <a^c< — ab. 

(a - b)(b - c) {a- b) 2 

On the other hand, using that < c < b < a, we get that 

2 111 111111 
c < A\ <b <=> < — = 1 < 

bacA2A\bacba 

2 111 b 11 

b a c b (a + c)b - ac b a 
ac 

<=> b > — . 

a + c - yfac 

Thus,ii e (0,c)and/l2 € (b, a) if and only if c < (a-2b)ab/(a-b) 2 andb > ac/(a+c- yjac). □ 
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Let us look for the winding numbers of the trajectories described in the previous proposition. 
The winding numbers mo, m\, and mi describe how the periodic billiard trajectories fold in K 3 . 



The following results can be found in [24, table 1]. First, mo is the period. Second, m\ is the 
number of xy-crossings and mi is twice the number of turns around the z-axis for EHl-caustics; 
m\ is twice the number of turns around the x-axis and mi is the number of yz-crossings for 
EH2-caustics; m\ is the number of tangential touches with each hyperboloid of one sheet caustic 
and mi is twice the number of turns around the z-axis for HlHl-caustics; m\ is the number 
of xz-crossings and mi is the number of yz-crossings for HlH2-caustics. Besides, all periodic 
trajectories with HlHl-caustics or HlH2-caustics have even period. Several periodic billiard 
trajectories with elliptic period m = 3 were depicted in j25l Table XV and Table XVII]. We 
conclude by direct inspection of those pictures that the nonsingular periodic billiard trajectories 
inside a triaxial ellipsoid with elliptic period m = 3 have winding numbers 

mi = 2, mi — 4, mo = 6. 

This agrees with the formulas (fTTb given in TheoremQj] We emphasize that those formulas were 
not rigorously proved, because their "proof" was based on ConjectureQ] 

Next, we establish the algebraic formulas for the caustic parameters of other nonsingular 
periodic billiard trajectories. We begin with a technical lemma about four-degree polynomials. 

Lemma 19. Let Q(x) = (x - a-)(x - fl-)(x - B + )(x - a+) for some a_ < B- < B + < a+. Let 
v_ e (ar_,j8_), v e (J3-,B + ), and v + e (B + , a + ) be the three roots of Q'(x). Then 

G(v+) < Q(v-) »ff- + ff + >iS- 

Proof If we set 77 = (J3+ + /3_)/2 and £ = {B+ -j3-)/2, then 

Qfo + s)- Q(j] - s) = 2(s 2 - f)(J3- +p + -a_- a + )s, \fs e M. 

On the one hand, if a_ + a+ > /L + then Q(j] + s) < Q(r\ - s) for all s > if, which implies that 
G(v+) < Q(v )■ On the other hand, if a_ + a+ < /3- + /3 + , then Q(j] + s) > Q(j] - s) for all s > 
which implies that Q{v + ) > Q{vJ). Finally, if a_ + a+ = >8_ + P+, then Q(t] + s) - Qirj - s) for 
all s e E, which implies that Q(v + ) = Q(v_). □ 

We can now answer some questions about the nonsingular periodic billiard trajectories found 
in the second item of Theorem[8] although the study is restricted to the spatial case. 

Proposition 20. There exist periodic billiard trajectories inside the triaxial ellipsoid M6$ with 
elliptic period m — 4, signature t — (0, 0, 1), and caustic type HI HI if and only if 

c < ab/(a + b). 

Besides, the caustic parameters A\ and Ai of such periodic billiard trajectories are the roots of 
the quadratic polynomial (s\ — si)t 2 /2 - sit + 1, where 

s, = l/a l + l/b l + l/c'-2/d', 1=1,2, (18) 

and d is the only root of the cubic polynomial t 1 - 2(a + b + c)t 2 + 3(ab + ac + bc)t - Aabc in the 
interval (a, +00). 
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Proof. Ifr = (0,0, 1), the decompositions presented in Definition|6]are 7 T = {2,3}, K T = {1,4,5}, 
V T = {1}, and W T = 0. Thus, C(4, 3; r) holds if and only if there exists some 5\ e (0, 75) such that 
the following two equivalent properties hold: 

(i) P(x) = {x- 6i) 2 (x - y 2 )(x - y 3 ) => Q(x) := x(x - 7l )(x - y A )(x - 75) = P(x) - P(0). 

(ii) y' 2 + y\ + 2S\ = y\ +y\ + y< 5 , for / = 1, 2, 3. 

If the caustic type is H1H1, then < c < Ai < A% < b < a, so 

71 = 1/c, 72 = lMi. 73 = IM2, 74 = 75 = 1/fl- 

Let e; = 6/(71,74,75) for / = 1,2,3. We set d = l/8\ > a. Then Q(x) = x 4 - e-ix 3 + q^x 1 - &\x 
and 6\ is a root of Q'(x) - Ax 1 - 3eijc 2 + 2e2^c - 63. Hence, d is a root of the cubic polynomial 

q(t) = -abet 3 Q '(1 /?) = f 3 - 2(a + + c)t 2 + 3(ab + ac + bc)t - Aabc. 

We note that g(0) = -Aabc < 0, q(b) = -b(b - a)(b - c) > 0, q(a) = -a(a - b)(a - c) < 0, and 
lim f ^ +co q(t) = +00. This shows that q{t) has just one root in the interval (a, +00). 

From property ([n} above, we deduce that the sums si := l/A\ + l/A^ verify relations ( fT8b . 
Besides, 1/A\ and I/A2 are the roots of (x - l/Ai)(x - I//I2) = jc 2 - siJc + (s 2 - «2)/2, so that A\ 
and /I2 are the roots of the quadratic polynomial (s 2 - S2)t 2 /2 — sit + 1, 

We look for ellipsoidal parameters such that the previous periodic trajectories exist. From 
property (0} above, we deduce that such ellipsoidal parameters exist if and only the graph {y - 
Q(x)} intersects the horizontal line {y = Q{5\)} at two different points 72, 73 e (74, 75) or, equiv- 
alent^, if and only if Q{Sj) < Q{6\), where S\ < 62 < 63 are the three ordered roots of the 
derivative of the polynomial Q{x) = x(x - y\)(x - ji)(x - 75). But 

Q(<%) < Q(8i) <=> 71 > 74 + 75 <=> c < ab/(a + b), 

according to Lemma[T9l □ 

As we have explained before, the period and winding numbers of any nonsingular periodic 
billiard trajectory can be determined by direct inspection of its corresponding figure. A periodic 
billiard trajectory with elliptic period m - A and caustic type H1H1 whose caustic parameters 
verify the relations given in Proposition l20l is displayed in [24, Figure 13]. That trajectory has 
(Cartesian) period mo = A and winding numbers 

m2 = 2, ni\ — 3, mo = 4. 

Hence, gcd(TOo, m\, mj) = 1, so the elliptic winding numbers are TO2 = 2,rh\ = 3, and too = 4. 
This result reinforces Conjecture[3] Besides, these nonsingular four-periodic billiard trajectories 
with caustic type H1H1 are quite interesting, because they display the minimal period among all 
nonsingular periodic billiard trajectories; see 11241 Theorem 1]. 

We end the study at this point. We just mention that there exist similar results when the 
signature or the caustic type do not coincide with the ones given in Proposition|20] Analogously, 
the case to = 5 can be dealt with using the same techniques, although the final formulas become 
more complicated. For instance, it can be easily checked that the caustic parameters A\ and 
A2 of the billiard trajectories with elliptic period to = 5 and signature r = (1,1,0) verify the 
homogeneous symmetric polynomial equations 

8s 3 + s] - 6siS 2 , I6S4 + s 4 - As]s 2 + As\, 
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where 5, = l/a l +l/b l +l/c l +l/A^+l/A[ for/ = 1,2,3,4. Each of the other signatures t € T(5, 3) 
gives rise to similar homogeneous — although not symmetric — polynomial equations of degrees 
three and four in the variables I /a, l/b, 1/c, l/Ai and 1/A2- We left the details to the reader. 
Finally, we remember that the original matrix formulation of the generalized Cayley condition 
C(5, 3) gives rise to two homogeneous symmetric polynomial equations of degrees 23 and 24 in 
those five variables, as explained in Section [3] This confirms, once again, that the polynomial 
formulation offers great computational advantages over the matrix formulation. 
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